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SUMMARY 


SYMBOLS 


The  theory  of  Lagrangian  multipliers  is  applied  to  ike  prob¬ 
lem  of  finding  both  upper  and  lower  limits  to  the  true  compressive 
buckling  stress  of  a  clamped  rectangular  plate.  The  upper  and 
lower  limits  thus  bracket  the  true  stress,  which  cannot  be  exactly 
found  by  the  differential-equation  approach.  The  procedure  for 
obtaining  the  upper  limit,  which  is  believed  to  be  new,  presents 
certain  advantages  over  the  classical  Rayleigh-Ritz  method  of 
finding  upper  limits.  The  theory  of  the  lower-limit  procedure  has 
been  given  by  Trefftz  but,  in  the  present  application,  the  method 
differs  from  that  of  Trefftz  in  a  way  that  makes  it  inherently  more 
quickly  convergent.  It  is  expected  that  in  other  buckling  prob¬ 
lems  and  in  some  vibration  problems  the  Lagrangian  multiplier 
method  of  finding  upper  and  lower  limits  may  be  advantageously 
applied  to  the  calculation  of  buckling  stresses  and  natural 
frequencies. 

INTRODUCTION 

Many  important  problems  that  cannot  be  exactly  solved  by 
the  differential-equation  approach  and  must  therefore  be 
analyzed  by  approximate  methods  arise  in  the  buckling  and 
vibrations  of  thin  plates.  The  theory  of  Lagrangian  multi¬ 
pliers  can  be  a  powerful  tool  in  the  analysis  of  many  of  these 
problems.  The  present  paper  presents  the  details  of  applica¬ 
tion  as  well  as  the  fundamental  principles  of  the  Lagrangian 
multiplier  method  by  demonstrating  the  use  of  the  method 
to  obtain  both  upper  and  lower  limits  to  the  true  compressive 
buckling  stress  of  a  rectangular  plate  clamped  along  all  edges. 

The  procedure  for  obtaining  the  lower  limit  is  similar  to  a 
method  used  by  Trefftz  (reference  1)  and  recently  described 
by  Reissner  (reference  2).  The  present  lower-limit  method 
differs  from  that  of  Trefftz,  however,  in  a  way  that  makes  it 
inherently  more  quickly  convergent.  The  upper-limit  pro¬ 
cedure,  which  does  not  appear  to  have  been  presented  pre¬ 
viously,  is  simpler  than  the  usual  Rayleigh-Ritz  method  and 
may  be  expected  to  permit  the  computation  of  more  accurate 
results  with  less  labor. 

In  a  recent  treatment  of  the  problem  of  compressive  buck¬ 
ling  of  clamped  plates,  extensive  calculations  of  lower  limits 
were  made  by  Levy  (reference  3)  by  means  of  a  procedure 
equivalent  to  the  Trefftz  method.  The  results  were  esti¬ 
mated  by  Levy  to  be  within  0.1  percent  of  the  true  results. 
In  order  to  illustrate  the  methods  of  the  present  paper,  upper 
and  lower  limits  to  the  buckling  stress  of  a  square  plate  are 
computed  to  within  0.1  percent  of  each  other;  a  positive 
check  on  the  accuracy  of  Levy's  results  is  thus  obtained. 


a  length  of  plate,  in  direction  of  stress 

b  width  of  plate,  perpendicular  to  stress 

/S  aspect  ratio  (a/b) 

t  thickness 

y  Poisson's  ratio 

E  Young’s  modulus  of  elasticity 


D  plate  stiffness  in  bending  (fizffZTffj) 


x  plate  coordinate  in  direction  of  stress 

y  plate  coordinate,  perpendicular  to  direction  of 

stress 

w  plate  buckling  deformation,  normal  to  plane  of 

the  plate 

<rs  critical  compressive  stress,  in  ^-direction 

k  critical  compressive  stress  coefficient  in  the 

formula  <rx=k 

V  internal  energy  of  deformation 

T  external  work  of  applied  stress 


an  Fourier  coefficient  of  cos  — - 

*  a 

br  Fourier  coefficient  of  cos 7-^ 

a 


am .  Fourier  coefficient  of  cos  — -  cos  — 

a  b 

i,j,m,n,p,q  even  integers 
r,  s  odd  integers 

Smn  Kronecker  delta  (1  if  m=n;  0  if  m^ri) 

A _ I _ 

[(m*+nsFy-km*l?]  (1+5*0+ 50n) 

a,  X,  \},  Ait,  -y  Lagrangian  multipliers 

THEORETICAL  BACKGROUND 

Rayleigh-Ritz  method. — The  Rayleigh-Ritz  energy  method 
for  determining  the  critical  stress  of  a  thin  plate  consists  of 
the  following  steps: 

(1)  The  deflection  surface  of  the  buckled  plate  is  expressed 
in  expanded  form  as  the  sum  of  an  infinite  set  of  functions 
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having  undetermined  coefficients.  In  general,  each  term  of 
the  expansion  must  satisfy  the  geometrical  boundary  condi¬ 
tions  of  the  problem. 

(2)  The  energy  of  the  load-plate  system  is  computed  for 
this  deflection  surface  and  is  then  minimized  with  respect  to 
the  undetermined  coefficients. 

(3)  This  minimizing  procedure  leads  to  a  set  of  linear 
homogeneous  equations  in  the  undetermined  coefficients. 
These  equations  have  nonvanishing  solutions  only  if  the 
determinant  of  their  coefficients  vanishes.  The  vanishing  of 
this  stability  determinant  provides  the  equation  that  may  be 
solved  for  the  buckling  stress. 

When  the  set  of  functions  used  is  a  complete  set  capable  of 
representing  the  deflection,  slope,  and  curvature  of  any  possi¬ 
ble  plate  deformation,  the  solution  obtained  is,  in  principle, 
exact.  Since,  however,  the  exact  stability  determinant  is 
usually  infinite,  a  finite  determinant  yielding  approximate 
results  is  used  instead. 

Lagrangian  multiplier  method. — The  Lagrangian  mul¬ 
tiplier  method  follows  the  general  procedure  outlined  for 
the  Rayleigh-Ritz  method,  with  but  one  outstanding  change. 
The  restriction  in  step  (1)  that  the  boundary  conditions  be 
satisfied  by  every  term  of  the  expansion  is  discarded  and  is 
replaced  by  the  condition  that  the  expansion  as  a  whole  sat¬ 
isfy  the  boundary  conditions.  This  condition  is  mathemat¬ 
ically  satisfied  in  step  (2),  during  the  minimization  process, 
by  the  use  of  Lagrangian  multipliers. 

The  fundamental  advantage  of  the  Lagrangian  multiplier 
method  lies  in  the  fact  that,  with  the  rejection  of  the  neces¬ 
sity  of  the  fulfillment  of  boundary  conditions  term  by  term, 
the  choice  of  an  expansion  is  much  less  restricted.  In  the 
clamped-plate  compression  problem,  a  simple  Fourier  expan¬ 
sion  may  be  used  instead  of  the  complicated  functions 
assumed  in  the  Rayleigh-Ritz  analyses  of  this  problem  (refer¬ 
ences  4  and  5).  Furthermore,  the  orthogonality  properties 
of  the  simple  Fourier  expansion  lead  to  energy  expressions 
of  a  simplicity  that  is  instrumental  in  permitting  accurate 
computations. 

Approximate  solutions:  upper  and  lower  limits, — 
The  Lagrangian  multiplier  method,  as  well  as  the  Rayleigh- 
Ritz  method,  gives  a  theoretically  exact  solution  for  the 
buckling  stress;  but  ordinarily  only  approximate  results  are 
obtained  because,  of  the  practical  necessity  of  considering 
finite  rather  than  infinite  determinants.  In  the  Rayleigh- 
Ritz  method  the  approximate  result  is  always  higher  than 
the  true  buckling  stress.  In  the  Lagrangian  multiplier 
method,  however,  it  is  possible  to  obtain  approximate  solu¬ 
tions  in  two  different  ways,  winch  permit  the  computation 
of  a  lower  limit  as  well  as  an  upper  limit  to  the  true  buckling 
stress.  As  determinants  of  higher  order  are  used  to  obtain 
approximations  of  higher  order,  both  the  upper-limit  and 
lower-limit  results  approach  the  true  buckling  stress.  Thus, 
the  Lagrangian  multiplier  method  can  be  used  to  provide  a 
result  to  within  any  specified  degree  of  accuracy.  It  may 
be  expected,  furthermore,  that  a  particular  determinant  in 
the  Lagrangian  multiplier  method  ought  to  yield  a  more 
accurate  result  than  a  determinant  of  equal  order  in  the 
Rayleigh-Ritz  method. 


LAGRANGIAN  MULTIPLIERS 

The  procedure  used  in  applying  the  fundamental  mathe¬ 
matical  principles  of  Lagrangian  multipliers  is  described  in 
this  section;  a  general  proof  of  the  validity  of  the  method  is 
given  in  the  appendix. 

Let  it  be  required  to  minimize  a  function  of  N  variables 

f(.XUXt,  X3,  .  .  .  xN)  (1) 

where  the  x’s  are  not  independent  but  are  bound  together 
by  the  relationship 

p(z i,  *2,  *3,  •  •  •  aw)=0  (2) 

Lagrange’s  method  of  simultaneously  minimizing  /  and 
satisfying  the  constraining  relationship  (2)  is  to  minimize 
the  function 


/-’-Xp 


with  respect  to  the  x’s.  The  quantity  X  is  the  undetermined 
Lagrangian  multiplier.  The  necessary  conditions  for  min¬ 
imizing/  then  become 


\g.-x<|_0  (K.1,2  A...JV) 

p= o  (2) 

Note  that  these  expressions  are  N+ 1  equations  in  the  N-\- 1 
unknowns  xu  Xj,  ...  xN  and  X. 

If  there  are  two  relationships  that  constrain  the  x’s; 
that  is,  if 

.  _  <Pi{Xi,X2,  ■  ■  •  zw)=0 


Pi(_X i,  X%t  .  .  .  Xjv)  —  0 


two  Lagrangian  multipliers  are  then  needed.  The  function 
to  be  minimized  becomes 


/ — XjySi — Xj^Oj 


and  the  minimizing  equations  are 


dyi 

dxjz 


-x,|fc- 

dxtr 


(#=1,2,3,  ...  AO 


<Pi=0 


P2=0 


The  method  is  easily  extended  to  cover  the  case  of  any 
number  of  constraining  relationships. 

PRELIMINARY  ILLUSTRATIVE  EXAMPLE 

Before  the  main  example  is  given,  a  simpler  buckling 
problem  will  be  analyzed  by  the  Lagrangian  multiplier 
method  in  order  that  the  method  of  application  of  Lagrangian 
multipliers  may  be  most  clearly  presented  without  the  ob¬ 
scuring  details  of  analysis  of  more  complicated  problems. 
This  elementary  problem  requires  the  use  of  but  a  single 
Lagrangian  multiplier,  which  leads  to  a  single  stability 
equation. 
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Consider  a  square  plate,  clamped  along  two  opposite 
edges,  simply  supported  along  the  other  two  edges,  which  is 
loaded  in  compression  on  the  simply  supported  edges.  (See 
sketch.)  From  the  exact  solution  of  this  problem,  the  deflec¬ 
tion  surface  of  the  plate  is  known  to  be  sinusoidal  in  the 
x-direction.  The  deflection  in  the  y-direction,  known  to  be 
symmetrical,  must  satisfy  the  elamped-edge  boundary 
conditions;  that  is,  zero  deflection 

w(x,0)=w(x,a)=0  (3) 

and  zero  slope 

||(x,0  )-£(«,  ®)=0  (4) 

The  present  method  uses  a  cosine-series  expansion,  whereas 
the  Trefftz  procedure  would  use  a  sine-series  expansion.  The 
problem  is  solved  by  both  methods  for  comparison. 


Cosine-series  solution. — In  the  cosine-series  solution  the 
expansion 

m 

.  mirx  VV  nry 

w=sm—  >  a„  cos  (5) 

n-oi2,4 

may  represent  the  deflection  surface  having  m  half-waves 
in  the  ^-direction,  since  the  Fourier  series  of  even  cosines  is 
a  complete  symmetrical  set. 

The  boundary  conditions  (equation  (4))  on  the  slope  are 
satisfied  by  each  term  of  the  expansion;  however,  in  order 
that  w  satisfy  the  conditions  of  equation  (3)  on  the  edge 
deflection,  it  is  necessary  that 

00 

Y  a*=0  (6) 

n^4 

Equation  (6)  is  a  constraining  relationship  on  the  a’s  and 
as  such  will  be  introduced  in  Lagrange's  minimization  process. 

As  in  the  Rayleigh-Ritz  method,  the  internal  energy  of 
deformation  and  the  external  work  of  the  stresses  are  then 
calculated.  Using  the  value  for  w  as  given  by  equation  (5) 


in  the  general  formulas  (reference  6,  equations  (199)  and 
(201)  (modified)) 


gives 


-rr_D  ra  f*  (  /d2w,tfw\2 
V-2j0  Jo 


(8) 


n=0,2.4 

The  usual  Rayleigh-Ritz  procedure  requires  that  the 
expression 

V—T  (9) 


be  a  minimum  with  respect  to  the  a’s.  In  the  present  ex¬ 
ample,  however,  the  a’s  are  not  independent  but  are  bound 
by  equation  (6).  Hence,  mathematically  stated,  the  ex¬ 
pression  V — T  must  be  a  minimum,  subject  to  the  constraint 
relationship  on  the  a’s 

oo 

Yj  a»=°  (6) 

nU0A4 

Solving  this  minimization  problem  by  Lagrange’s  method 
makes  it  necessary  to  minimize 

(10) 

»-0A4 


with  respect  to  the  a’s.  The  necessary  conditions  for  a 
minimum  then  become 


b{V-T) 


(J-  0,2,4,...) 


*=0,2,4 


a»=0 


or,  upon  differentiation  and  simplification, 

4/7^ 

Cl+5or)[(m2+ji),-m*i]a,-^)X=0  (j=0,2,4,  .  .  .)  01 


Y  a»=0  (0) 

n-0A4 
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Solving  equation  (11)  for  as  and  substituting  into  equa¬ 
tion  (6)  gives  the  stability  equation  that  determines  k: 


Yj  [(m2+j2)s— m^Kl+^r0  (12) 

Z-0,2,4 

For  a  particular  number  of  half  waves  m,  this  equation 
may  be  solved  by  evaluating  the  series  for  several  trial 
values  of  k  and  interpolating  to  find  the  k  that  makes  the 
series  vanish.  The  correct  value  of  m  is  that  which  gives 
the  lowest  value  of  k.  For  two  half  waves  (m= 2)  in  the 
loaded  direction,  the  theoretically  exact  value  of  k= 7.69 
(reference  6,  p.  345)  is  obtained  when  only  10  terms  of 
equation  (12)  are  computed. 

Sine-series  solution  (Trefftz  method). — The  same  problem 
will  be  treated  in  the  manner  suggested  by  Reissner  (refer¬ 
ence  2),  which  is  similar  to  Trefftz's  method  (reference  1). 

Let 

a> 

.  mvx  ,  .  rry 

w=sm  —  2j  hr  sm  a  (*3) 

The  boundary  conditions  on  deflection  (equation  (3))  are 
now  satisfied  term  by  term,  but  the  conditions  on  the  edge 
slopes  (equation  (4))  are  satisfied  only  by  making 


(14) 


Now,  the  expression  V—  T  is  computed  from  formulas  (7) 
and  (8)  by  using  the  value  of  w  given  by  equation  (13); 
then,  by  application  of  Lagrange’s  procedure, 


ao 

( V-T)-y ^  rbr 


(15) 


r-1,8.5 


must  be  minimized  with  respect  to  b,  where  s=l,  3, 5,  .  .  .. 
The  Lagrangian  multiplier  is  y. 

The  minimization  equations  are 


stability  equation,  the  Lagrangian  multiplier  method  is 
preferably  used  to  satisfy  the  zero-deflection  condition  rather 
than  the  zero-slope  condition.  Slope  is  the  derivative  of 
deflection,  and,  in  general,  differentiation  of  a  Fourier  series 
makes  it  more  slowly  convergent. 

THE  COMPRESSIVE  BUCKLING  OF  A  RECTANGULAR  PLATE 
CLAMPED  ALONG  ALL  EDGES 

The  previous  elementary  example  required  only  a  simple 
Fourier  "expansion  and  but~one  Lagrangian  multiplier  to  sat¬ 
isfy  the  boundary  conditions.  The  more  difficult  problem 
of  finding  the  buckling  stress  of  the  rectangular  plate 
clamped  on  all  edges  and  loaded  as  shown  in  the  accompany¬ 
ing  sketch  necessitates  a  double  Fourier  series,  as  well  as  an 
infinite  set  of  Lagrangian  multipliers  to  satisfy  the  boundary 
conditions. 


Boundary  conditions. — The  boundary  conditions  of  the 
problem  are: 

Zero  deflection,  loaded  edges 

w(0,  y)=v>(a,  y)= 0  (18) 


[(m2-!-®2)2— m2k]b,— 


4a2 


ys=0  (s=  1, 3, 5,  .  .  .)  (16) 


2  (14) 

Solving  equation  (16)  for  b,  and  substituting  in  equa¬ 
tion  (14)  gives  as  the  stability  equation 


co 

£ 

*-1,8,8 


(m2+82)2—m2k 


(17) 


Comparison  and  discussion  of  results. — The  series 
in  equation  (17)  converges  approximately  as  1/a2,  whereas 
the  series  in  equation  (12)  converges  approximately  as  1  fj*. 
Because  of  the  more  rapid  convergence  obtained  in  the 


Zero  deflection,  unloaded  edges 

w(x,  0)=w(x,  b)—0  (19) 


Zero  slope,  loaded  edges 

§>,lrt=f£(a,y>=0  (20) 

Zero  slope,  unloaded  edges 

^  (x,  0)=^  (x,  b)  =  0  (21) 

Fourier  expansions. — In  order  to  achieve  a  rapidly 
convergent  solution,  the  principles  established  by  the  pre¬ 
ceding  example  are  used  as  the  basis  for  choosing  the  Fourier 
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expansions  to  satisfy,  term  by  term,  the  conditions  of  zero 
slope  rather  than  those  of  zero  deflection. 

The  buckling  deformation  corresponding  to  the  lowest 
buckling  stress  is  always  symmetrical  perpendicular  to  the 
direction  of  load  but,  depending  on  the  aspect  ratio  of  the 
plate,  may  be  symmetrical  or  antisym  m  etrical  in  the  direc¬ 
tion  of  load.  Thus,  for  symmetrical  buckling,  let 


Since  this  Fourier  series  must  vanish,  each  infinite  series  that 
constitutes  a  coefficient  of  a  cosine  term  must  vanish.  (All 
the  Fourier  coefficients  of  the  Fourier  expansion  of  the  func¬ 
tion  zero  are  zero.)  Hence 

oo 

^  ani=  0  C?=0, 2, 4,  .  .  .)  (25) 

ra=0, 2,4 


cd  m 

EV\  mxx  nxy 

>  aan  cos  —  cos-g^  (22) 

ra— 0,2, 4  n= 0,2,4 

and,  for  antisymmetrical  buckling,  let 

w=  Y  cos  ^ cos  (23) 

rirte 

It  is  sufficient,  for  purposes  of  demonstration,  to  consider 
only  the  case  of  symmetrical  buckling.  Hereinafter,  w  there¬ 
fore  refers  to  the  value  given  by  equation  (22). 

Energy  expressions. — Using  the  expansion  given  by 
equation  (22)  in  the  evaluation  of  the  general  energy  and 
work  integrals  of  equations  (7)  and  (8)  gives 

<»  oo 

^  ^  (mHnW(l-M*o+S0„)<W2 

m-0,2,4  n-0,2,4 


By  expressing  the  fact  that  there  is  zero  deflection  along 
the  unloaded  edges  (equation  (19)),  it  can  be  similarly  shown 
that 

oo 

2  O'tn—O  (i=0, 2,4,  .  .  .)  (26) 

Now,  V—  T  must  be  a  minimum  with  respect  to  the  a’s, 
which  are  bound  by  equations  (25)  and  (26).  As  the  prob¬ 
lem  now  stands,  however,  it  is  not  in  the  form  to  which 
Lagrange’s  minimization  process  can  be  applied  since  V—  T_ 
does  not  contain  ctoo,  whereas  the  constraint  relationships 
do  contain  a^.  Hence,  Ooo  is  eliminated  from  the  con¬ 
straint  relationships  by  subtracting  the  first  of  equations 
(26),  the  equation  for  i=0,  from  the  first  of  equations  (25), 
the  equation  for  j— 0.  The  final  set  of  necessary  con¬ 
straining  relationships  on  the  minimization  of  the  energy 
expression  (24) 


T 


CD  CO 

v\  V 

8 a  Zj  Zj 

n=0, 2, 4  n=0,  2,4 


ml(l+5o»)a»»* 


Then 

v-t^-3 

8a® 

where 

~t— — [ (m*-{-  n’/S2)2— ( 1+ 5„o + 5o»)  ■ 
Note  that  V—  T  is  independent  of  am,  since  -^=0. 


Y  Y  Am*a’ 


(24) 


m-0,2,4  n-0,2,4 


Constraining  relationships. — The  boundary  conditions 
of  zero  slope  (equations  (20)  and  (21))  are  satisfied  by 
each  term  of  the  expansion  of  equation  (22),  but  the  condi¬ 
tions  on  deflection  (equations  (18)  and  (19))  must  be  satisfied 
by  the  expansion  as  a  whole.  Subs  titutingu? in  to  equation  (18) 
gives,  along  the  loaded  edges, 


18=0,2,4 


OO  CO 

W(0,  2/)=«?(u,  y)='^J  OaO+COS^^  a  m2 

1=0,2, ' 

CD 

E 

b-0,2 

+  COS^2  Cm, 


ra- 0,2,4 

-fcos^  2_J 

ra-0T2,4 


+ 


n-0,2,4 

=0 


then  becomes 


co 

QmO —  (Zo#=0 

=2^4,6 

03 

^  am}=0 

.^072,4 

CO 

E 


to -2, 4,6  n=2,4,6 


(i= 2,4,6,...) 


0  (*=2,4,6,...) 


n-0,2,4 


y  (27) 


Theory  of  upper  and  lower  limit  solutions. — A  theoretically 
exact  solution  to  the  problem  would  be  obtained  if  the 
energy  expression  (equation  (24))  were  minimized  with 
respect  to  all  the  a’s  and  at  the  same  time  all  the  rela¬ 
tionships  (27)  were  satisfied.  This  result  follows  from  the 
facts  that:  (a)  the  expansion  of  equation  (22)  is  a  complete 
symmetrical  set,  capable  of  representing  the  exact  symmet¬ 
rical  buckling  deformation,  and  (b)  the  fulfillment  of  the 
conditions  of  equations  (27)  ensures  that  the  boundary  con¬ 
ditions  are  completely  satisfied.  An  exact  solution  is  not 
possible,  however,  because  it  would  involve  an  infinite  deter¬ 
minant,  so  that  two  different  modifications  of  the  ideal  pro¬ 
cedure  are  used  to  obtain  approximate  results.  One  of  these 
methods  gives  an  upper  limit  to  the  true  buckling  stress, 
whereas  the  other  gives  a  lower  limit. 
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An  upper  limit  to  the  buckling  stress  can  be  found  by 
arbitrarily  setting  some  a’s  equal  to  zero,  minimizing  expres¬ 
sion  (24)  with  respect  to  the  remaining  a's,  and  satisfying  all 
the  constraint  relationships  (27) .  An  upper  limit  is  obtained 
inasmuch  as  arbitrarily  setting  some"  of  the  Fourier  co¬ 
efficients  equal  to  zero  has  the  effect  of  restraining  the  deflec¬ 
tion  of  the  plate,  which  in  effect  stiffens  the  interior  of  the 
plate  and  increases  the  stress  required  to  buckle  it. 

A  lower  limit  to  the  buckling  stress  can  be  found  by 
minimizing  expression  (24)  with  respect  to  all  the  a’s  but 
satisfying  only  some  of  the  constraining  relationships  (27). 
Neglecting  some  of  the  constraining  relationships  has  the 
effect  of  giving  the  plate  greater  freedom  at  the  edges  and 
hence  reducing  the  stress  required  to  buckle  the  plate. 

Lower  limit  solution. — In  accordance  with  the  require¬ 
ments  for  a  lower  limit,  the  constraining  relationships  (27) 
will  be  satisfied  only  up  to  j=q  and  %—p.  By  Lagrange’s 
minimization  process,  the  function  to  be  minimized  is  then 


6- 


ir^Db 

a~83r 


OP  00  /  CO  CO  \ 

X  2  2  a°”) 

i-G.2,4  n-0,2,4  \m-2,4,Q  n-2,4,6  / 


n  oo  p  co 

“  Zj  X'  2  am,~  Li  2  ata 


1- 2,4,6  m=0,2,4 


f-2,4,0  n-0,2,4 


The  a,  X2  ■  .  .  X„,  Ms  .  •  •  are  Lagrangian  multipliers. 
The  equations  for  minimizing  V—  T  with  the  constraining 
relationships  (27)  on  the  a’s  then  become 


bG 

bamM 


0 


(m,n= 0,  2,  4, 


Equations  (27)  up  to  i=p,j=q 


(28) 


Evaluation  of  bGjbamn  gives 

“(So*— 5«io) — X»— Pm—0  (29) 

where  X„  appears  only  if  2gn^g  and  y.m  appears  only  if 
2  £m£p.  From  equation  (29) 


a  aP+d, «+« — 0  (d,  a — 2,4,6,  .  .  .) 

fl-f* 

Hence  for  any  particular  d,  e,  either 


or 


0 


A^t,  *+« 

s+«=Q 


The  first  alternative,  however,  ordinarily  would  require  k  to 
be  very  high,  corresponding  to  the  budding  stress  of  a 
buckling  mode  with  many  waves  in  both  directions.  For  the 
lowest  buckling  load,  then, 


ctp-td.  !+.= 0  (d,e= 2,4,6,  .  .  .) 

It  is  therefore  necessary  to  be  concerned  with  only  the  other 
a’s,  which,  from  equation  (29),  are 


_  4 a*  .  ,  ,  . 

& mn  Pm) 

<lm0~TriDb  + 

4 

O' On  ^.4  JJ-fo 


(m,n?^o)  (30) 


In  equations  (30),  X„  does  not  appear  if  »>g  and  pm  does 
not  appear  if  m>p. 


Substituting  the  values  of  a’s  given  by  equations  (30)  back  into  the  constraining  relationships  (27)  up  to  j=g,  i=p  gives 

P 


\m-2,4,6  ft— 2, 4, 6  /  ro-2,4,6  n-2,4,6 

—  A<ija+  ^  ^  Amj=  0  (j=2,4,6,  .  .  .  q) 


w-2, 4, 6 


m-0.2,4 


oo  Q 

A(sXb=0  (1=5^4, 6, ...  p) 

n-lj,  2, 4  n=2,4.0 


(31) 


These  equations  form  a  set  of  Kp+2)  +  l  linear  homo¬ 
geneous  equations  in  a,  fa  .  .  .  p„,  X*  .  .  .  Xs.  Since  when 
buckling  occurs  the  a’s  are  not  all  zero,  by  equations  (30), 
the  Lagrangian  multipliers  are  not-  all  zero.  In  order  that 
equations  (31)  be  compatible,  the  determinant  of  the  coeffi¬ 
cients  of  the  Lagrangian  multipliers  must  vanish.  The 
vanishing  of  this  stability  determinant  provides  the  deter¬ 
minants!  equation  that  may  be  solved  for  k  by  substitution 
of  trial  values  and  interpolation. 


That  certain  elements  of  the  determinant  consist  of  an 
infinite  series  of  Amn  terms  is  evident;  these  series  convex’ge 
rapidly.  Since  such  rapidly  convergent  series  are  cal¬ 
culable  txrany  degree  of  accuracj7,  they  may  be  considered 
as  known  quantities.  Each  value  of  Amn  represents  the 
potential-energy  contribution  of  a  term  in  tho  expansion  for 
w;  hence,  the  effects  of  infinite  subsets  of  expansion  terms 
enter  into  this  solution.  Thus,  for  p=q= 2,  the  expansion 
terms  corresponding  to  the  a’s  shown  in  table  I  enter  into 
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the  solution;  similarly,  the  terms  represented  in  table  II 
enter  into  the  solution  when  p=q=4. 


TABLE  I.— FOUR  INFINITE  STRIPS  OF  FOURIER 
COEFFICIENTS  OF  EXPANSION  TERMS 


N'sv  m 

Tl  \ 

0 

2 

4 

e 

8 

.  .  . 

0 

an 

aw 

aw 

2 

on 

an 

041 

a« 

an 

. 

4 

<t2i 

6 

a« 

an 

S 

OflS 

aw 

TABLE  II.— SIX  INFINITE  STRIPS  OF  FOURIER 
COEFFICIENTS  OF  EXPANSION  TERMS 


\  ZB 

Zl  \ 

• 

2 

4 

6 

8 

•  -  - 

0 

at« 

an 

041 

aw 

aw 

2 

003 

an 

041 

a« 

au 

. 

4 

aoi 

aj< 

044 

a« 

a»i 

.  .  . 

6 

m 

an 

04C 

8 

aot 

an 

04J 

Upper-limit  solution. — The  lower-limit  solution  satisfied 
only  some  of  the  constraining  relationships  (27)  but  assumed 
the  existence  of  all  the  Fourier  coefficients.  If  an  upper 
limit  is  to  be  obtained,  it  will  be  necessary  to  satisfy  all  the 
constraining  relationships  while  arbitrarily  assuming  some 
fl’s  to  be  zero. 

As  a  direct  result  of  the  necessity  of  satisfying  all  the 
constraining  relationships  in  the  upper-limit  solution,  it  is 
found  that  the  first  of  equations  (27)  is  redundant  and  may 
be  discarded,  since  it  is  automatically  satisfied  when  all  the 
remaining  equations  (27)  are  satisfied.  As  a  proof  of  this 
redundancy,  the  conditions 

co 

amj=  0  .  (7  =  2,4, 6,  .  .  .) 

ra-0.2,4 

are  summed  over  j  and  subtracted  from  the  sum  of  the 
conditions 

o» 

^  (i=2, 4, 6,  .  .  .) 

'  a-oTit 

over  i  and  give 


It  is  to  be  emphasized  that  the  redundancy  of  a  constrain¬ 
ing  relationship  is  a  peculiarity  of  only  the  upper-limit  solu¬ 
tion,  since,  as  shown  by  the  proof  given,  the  redundancy 
depends  on  the  fact  that  all  the  constraining  relationships 
must  be  satisfied. 

"With  the  elimination  of  the  redundant  condition,  the 
necessary  constraint  relationships  become 

CD 

^  amj=0  (7=  2, 4, 6,...)  (32) 

Ttt— 0^2, 4 

CD 

^«..= 0  (i=2, 4, 8,...)  (33) 

n-0T2.4 

At  this  point,  in  accordance  with  upper-limit  theory,  it 
is  necessary  arbitrarily  to  set  certain  a's  equal  to  zero.  It  is 
possible  to  take  advantage  of  the  Lagrangian  multiplier 
method  by  allowing  infinite  rather  than  finite  sets  of  a’s 
to  exist  and  still  to  obtain  a  stability  determinant  of  finite 
order.  Thus,  infinite  strips  of  coefficients  of  the  type  shown 
in  tables  1  and  2  can  enter  into  the  solution.  In  the  lower- 
limit  case,  the  existence  of  all  coefficients  was  assumed,  but 
the  coefficients  ap+d> 5+,  were  proved  to  be  zero;  in  this 
upper-limit  solution,  it  will  be  arbitrarily  assumed  that 
these  same  a’s  are  zero;  thus 

qj>+<i.  s+«=0  (^i  6=2,  4,  6,  ...  ) 

The  constraining  relationships  (32)  and  (33)  become 


(Imj — 0 

(7=2,4,  6, ...  2) 

(32a) 

amj= 0 

0=2+2,  g+4,  . .  .co) 

(32b) 

a(n=0 

(i=2,  4,  6,  ...  p) 

(33a) 

atn=  0 

(i=p+2,  p+4,  ...  co ) 

(33b) 

71=0,  2,  4 

The  function  to  be  minimized  is 


r  E  E  £  X/  E  ^ 

771=0,2,4  n=0,  2,4  /= 2,4,8  171-0,2,4  ■ 

CO  P  P  00 

■  E  x*  E  am}~  E  -  E  a 

/=j+2,  s+4  za— 0,2,4  f-2,4,  6  n=0, 2,  4 

v\  ^ 

2j  ^  2j  au 

(-p+2,p+4  71=0,  2,  4 


/“4+2,0+4  za— 0,2,4  f-2,4,  6  n=0, 2, 4 
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The  first  double  summation  of  this  equation  extends  over 
only  the  values  of  m  and  n  such  that 

m^p  if 


Setting 
to  exist, 


if  m>p 

0  then  gives,  for  all  the  a’ s  arbitrarily  allowed 


4 a8  rA 

^m[A” 


»(^nH“  A*  »i)  I 


where  Xo  and  m  do  not  exist.  Substituting  back  into  the 
constraint  equations  (32a),  (32b),  (33a),  and  (33b)  gives 


CO 

Ami+ 

CO 

:  0  (i=2, 4, 6,  .  . 

•  2) 

(34a) 

TO=0,  2,  4 

m~2,4,6 

P 

Amj- 1- 

-A-mlPm— 

1 0  (i=2+2,g+4,. 

•  •“) 

(34b) 

m=0,2,4 

m-2,4,0 

£ 

AjttXB+jw< 

Ain=° 

ccT 

<n 

ii 

•c* 

P) 

(35a) 

n-sue 

71=0,2,4 

£ 

£  A«= 0 

(i=p+2,p+4, .  . 

.00) 

(35b) 

n=2;4,0 

B-0,  X  4 

These  equations  involve  all  the  Lagrangian  multipliers. 
They  can  be  reduced  to  a  set  of  equations,  however,  in 
Xa  ...  Xs,  ^  ...  /I*  in  the  following  manner: 

From  equation  (35b),  for  i=p+2,  p+4,  .  .  . 


n-2.4,0 


n-0,2,4 


From  equation  (34b),  forj=g+2,  g+4,  .  .  . 


i  CT~2.4.8 


m=0,2,4 


Substituting  these  expressions  for  nt  and  X/  into  equations 
(34a)  and  (35a),  respectively,  gives  as  the  final  stability 
equations: 


nt  — 0,2,4  m=2,4,0 


=0  (i=2,  4,  6,  .  .  .  q) 


g  oo 

»= 2,4,6  i*-£F5«+4 


TO- 2,4,6 


£ 


■44m* 

V.  m— 0T2.4  »/ 


eo 

^  ^  -d.jn=0  (i=2,  4,  6, 


P) 


n-0,2,4 


(30) 


Equations  (36)  form  a  set  of  ^  (p+g)  linear  homogeneous 

equations  in  Xa  .  . .  Xe,  jus  . . .  ju,.  The  stability  determinant 
is  the  determinant  of  the  coefficients  of  the  X's  and  p's. 

It  is  of  interest  to  note  that  in  the  usual  Rayleigh-Ritz 
solutions  only  finite  sets  of  expansion  terms  are  ever  taken 
into  account,  and  the  order  of  the  determinant  obtained  is 
ordinarily  equal  to  the  number  of  terms  considered.  It  is 
then  reasonable  that  a  particular  determinant  obtained  by 
the  Lagrangian  multiplier  method,  which  considers  infinitely 
more  expansion  terms  than  a  Rayleigh-Ritz  determinant  of 
equal  order,  may  be  expected  to  give  a  more  accurate  result. 

Numerical  example. — For  the  case  of  a  square  plate,  (3=1, 
upper  and  lower  limits  were  computed.  The  results  for 
the  buckling-stress  coefficient  4  were:. 


Approximation 

Lower  limit 

Upper  limit 

9.  99 

10. 11 

10.08 

Second;  p=g=4. 

. - . 

— 

10.07 

The  expectation  that  the  Lagrangian  multiplier  method 
should  give  closer  upper  limits  than  the  Rayleigh-Ritz 
method,  for  a  given-order  determinant,  can  be  confirmed  for 
this  example.  A  second-order  Lagrangian  multiplier  deter¬ 
minant  gives  an  upper  limit  of  k=  10.11,  whereas  Maulbetscb 
(reference  4)  and  Smith  (reference  5)  use  complicated  de¬ 
flection  functions  in  the  Rayleigh-Ritz  method  to  derive 
third-order  determinants  that  give,  respectively,  4=10.48 
and  4=10.11. 

It  is  seen  that  the  second  approximation,  requiring  the 
evaluation  of  a  fourth-order  determinant  for  the  upper  limit 
and  a  fifth-order  determinant  for  the  lower  limit,  definitely 
establishes  the  value  of  4  to  within  0.1  percent. 

Levy  (reference  3)  used  an  ingenious  method  of  obtaining 
lower  limits  that  is,  in  fact,  equivalent  to  the  Trofftz  method 
of  using  double  sine  series  and  satisfying  the  zero  edge-slope 
condition  by  the  Lagrangian  multiplier  method.  On  the 
basis  of  computations  involving  determinants  up  to  order 
twenty,  Levy  concluded  that  his  results  obtained  from  tenth- 
order  determinants  are  within  0.1  percent  of  the  true  results. 
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inasmuch  as  Levy  obtained  £=10.074  for  the  square  plate, 
the  present  relatively  simple  upper-  and  lower-limit  calcula¬ 
tions  show  that  his  estimated  limit  of  error  is  correct  for 
this  case. 

CONCLUSIONS 

1.  The  Lagrangian  multiplier  method  can  be  used  to 
compute  accurate  upper  and  lower  limits  to  the  compressive 
bucklipg  stress  of  a  clamped  rectangular  plate,  thereby 
bracketing  the  true  result. 

2.  From  a  consideration  of  rapidity  of  convergence  toward 
the  exact  solution  in  clamped-plate  problems,  it  is  preferable 
to  use  an  expansion  that  satisfies  the  zero-slope  boundary 
conditions  term  by  term  rather  than  the  zero -deflection 
boundary  conditions. 

3.  Because  of  the  fact  that  the  Lagrangian  multiplier 
method  permits  the  effects  of  infinite  subsets  of  expansion 


terms  to  enter  into  the  solution,  it  is  believed  that  a  particular 
stability  determinant  derived  by  the  Lagrangian  multiplier 
method  will,  in  general,  yield  a  closer  upper  limit  than  that 
obtained  from  a  determinant  of  equal  order  in  the  Rayleigli- 
Ritz  method. 

4.  It  is  expected  that  the  method  of  Lagrangian  multipliers 
may  be  useful  in  the  analysis  of  other  stability  and  vibration 
problems.  In  particular,  the  method  may  be  immediately 
applied  to  the  determination  of  vibration  frequencies  of 
clamped  plates,  and  to  the  determination  of  buckling  stresses 
of  clamped  plates  under  compression  in  two  directions. 
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Langley  Field,  Va.,  May  S,  1946. 


APPENDIX 

GENERAL  PROOF  OF  THE  METHOD  OF  LA  GRAN GIAN  MULTIPLIERS 


Let  it  be  required  to  minimize 

f(xu  is,  xt,  ...  xN)  (Al) 

where  the  N  i’s  are  bound  by  the  P  independent  relation¬ 
ships  (P<iV) 

pj  (#1)  Xg,  X3,  . .  .  av)  =  0  («7=1,  2,  3 ,  ...  P)  (A2) 
It  will  be  proved  that  the  equations  for  determining  the 


minimizing  values  of  the  x's  are: 

Ji+x,H+XllS+'  •  ■+x'iS“°  CX-W.-JA  (A3) 

<Pj{xllxi,Xi,...xN)=0  (*7=1,  2,  3, .  .  P)  ■  (A2) 


The  X's  are  Lagrangian  multipliers;  these  JV+P  equations 
determine  N  i’s  and  P  X’s. 


If  the  values  of  only  N—P  x’s  are  known,  the  remaining  P  i’s  are  determined  from  the  P  relationships  (A2).  For 
convenience,  consider  the  last  P  x’s  in  equation  (Al)  to  be  dependent  upon  the  first  N—P  x's.  Then  for  /  to  be  a 
minimum  its  first  partial  derivatives  with  respect  to  the  independent  x's  must  vanish,  or: 


¥  .  £>/  daw-p+i  1  ¥  dZiy-p+2  ,  ¥_ 

bxM  axw_P+1  aijw  ai^_p+a  ai^  dxy  ax^ 


(M=l,  2,  3,  .  .  .  (N-P)) 


(A4) 


But  each  of  these  equations  contains  P  quantities  that  cannot  be  directly  evaluated — the  derivatives  of  the  depend¬ 
ent  variables  with  respect  to  the  independent  variables.  For  each  value  of  HZ,  these  P  derivatives  are  determined  by 
differentiating  each  of  the  P  constraint  relationships  (A2)  with  respect  to  xM.  Thus, 


a  <pj 


a ipj  ax^r-p+i ,  b<pj  axjv_p+2 1  i  dpj 

1  l  •  •  •  1 


?>Xm 


=0 


(J=  1,  2,  3,  .  .  .  P) 


(AS) 


1  dJV—P+I  1  P+2  QXjyf 

Now,  for  each  particular  value  of  M,  equation  (A4)  and  the  P  equations  (AS)  make  up  a  set  of  P+1  linear  homo¬ 
geneous  equations  in  the  P+1  quantities  1,— ^ 


da^p+i  bxN  P+i>  _  gjnce  tjiege  quantities  are  surely  not  all  zero, 


dx 


■M 


the  determinant  of  their  coefficients  must  vanish.  Hence,  it  is  found  that  for  j  to  be  a  minimum  it  must 
be  true  that: 


2L 

¥ 

¥  . 

..3L 

dxM 

dXy-p+i 

dijv-p+a 

dXy 

<V 

dpi 

dpi 

dpi 

dxM 

i)Xtr-p+i 

dXtr-p+2 

dxN 

dp* 

dpg 

d<p2 

= 0  (M=  1,2,3, 

dxM 

dxN_P+i 

dlAT-P+2 

dXy 

dpp 

dpp 

Si fp 

dpp 

dxM 

&Zy-p+ 2 

dxN 

necessarily 


(AG) 


It  will  now  be  demonstrated  that  these  necessary  minimization  equations  will  hold  if  eauations  (A3)  bold.  Interchanging 
the  rows  and  columns  of  the  determinant  in  equation  (A6)  gives: 


¥ 

dpi 

d<p2 

dlfip 

H* 

/O 

dxM 

dxu 

dxM 

¥ 

dipt 

d<pp 

dXy_p+i 

dXtf-P+l 

dXy-p+ 1 

dx^_P+1 

¥ 

dpi 

dlfii 

dipp 

dXy-p+a 

dXy_p^  2 

dXy-p+2 

dX.y_p+2 

¥ 

dpi 

d<f>2 

d<pp 

dXy 

bly 

bXy 

dx* 

(Af=I,2,  3, . . .  (A7— P)) 


(A  7) 
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The  vanishing  of  this  determinant  is,  however,  precisely  the  condition  of  compatibility  of  the  equations 

+  Xl^i  +  x,|^  +...+  =ol 

fa*  bxM  oxx 


_  y 


+Xi  - hXs 


^_+...+Xp?r^-=o 

oa^-p-fi  oa;^-p  vi 


sr^-+1''53re-+x=srs-+  ■  ■  •  +x's#2£—° 

OX.V-P42  O^iV-P+2  OZ&-P+2  0*2IT-P+ 2 


2)/  I  i  ,  ..  2  ,  ,  .  5^p  _n 

=r —  i  Ai  5: —  *r*  Ag  5: —  “r  •  •  •  ~r  Ap  — ~ —  —  U 

bxN  'bxt,  'b  xN  <pxK 


when  they  are  considered  as  linear  homogeneous  equations  in  the  quantities  1,  Xlf  X2,  .  .  .  \P. 

Since  a  determinant  (A7)  exists  for  each  value  of  M  up  to  (N — P),  a  set  of  equations  (A8)  exists  for  each  M.  It  is 
seen  that  in  these  sets  only  the- first  equation  varies,  since  only  the  first  equation  depends  upon  M.  Observation 
shows  that  all  the  ( N — P )  determinants  of  equation  (A7)  can  be  derived  from  the  set  of  N  equations 


e-+Xi: 


+  .  .  .+Xpl^=0  (i£=l,  2, 3, 


by  successively  writing  the  determinants  of  compatibility 
of  the  last  P  equations  with  each  of  the  first  ( N—P )  equa¬ 
tions  in  turn.  It  has  thus  been  proven  that  if  equations  (A3) 
are  true,  the  minimizing  equations  (A6),  equivalent  to 
equations  (A7),  must  hold. 

It  is  seen,  however,  that  equations  (A3)  are  N  equations 
in  ( N+P)  unknowns  consisting  of  N  x’s  and  P  X’s.  The 
remaining  necessary  P  equations  come  from  the  original 
equations  of  constraint  (A2).  Hence,  the  simultaneous 
equations  (A2)  and  (A3) 


5/  ! 2)^  ,  b<pj  . 
d^+Xl  dzK+ 


+Xp^=0  (#=  1,2,3, 


<pAx Xt,  Xt,  . . .  xh)=Q  (J=l,  2,  3,  ...  P)  (A2) 


are  necessary  equations  for  the  minimization  of  J(xu  x2,  xs, 
. . .  x^),  which  was  to  be  proved. 
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